Linear Differential Equations 
with Fuzzy Boundary Values 



Nizami Gasilov 

Baskent University, 

Eskisehir yolu 20. km, Baglica, 

06810 Ankara, Turkey 
Email: gasilov@baskent.edu.tr 



§ahin Emrah Amrahov 

Ankara University, 

Computer Engineering Department, 

06100 Ankara, Turkey 

Email: emrah@eng.ankara.edu.tr 



Afet Golayoglu FatuUayev 

Baskent University, 

Eskisehir yolu 20. km, Baglica, 

06810 Ankara, Turkey 

Email: afet@baskent.edu.tr 



o 



< 



> 



o 



X 



Abstract — In this study, we consider a linear differential 
equation with fuzzy boundary values. We express the solution 
of the problem in terms of a fuzzy set of crisp real functions. 
Each real function from the solution set satisfies differential 
equation, and its boundary values belong to intervals, determined 
by the corresponding fuzzy numbers. The least possibility among 
possibilities of boundary values in corresponding fuzzy sets is 
defined as the possibility of the real function in the fuzzy solution. 

In order to find the fuzzy solution we propose a method based 
on the properties of linear transformations. We show that, if 
the corresponding crisp problem has a unique solution then the 
fuzzy problem has unique solution too. We also prove that if the 
boundary values are triangular fuzzy numbers, then the value of 
the solution at any time is also a triangular fuzzy number. 

We find that the fuzzy solution determined by our method 
is the same as the one that is obtained from solution of crisp 
problem by the application of the extension principle. 

We present two examples describing the proposed method. 

Keywords: fuzzy boundary value problem, fuzzy set, linear 
transformation. 



I. Introduction 

Approaches to fuzzy boundary value problems can be of two 
types. The first approach assumes that, even if only the bound- 
ary values are fuzzy in the handling problem, the solution is 
fuzzy function, consequently, the derivative in the differential 
equation can be considered as a derivative of fuzzy function. 
This derivative can be Hukuhara derivative, or a derivative in 
generalized sense. Bede 1 1 J has demonstrated that a large class 
of boundary value problems have not a solution, if Hukuhara 
derivative is used. To overcome this difficulty, in f2l and fSl 
the concept of generalized derivative is developed and fuzzy 
differential equations have been investigated using this concept 
(see also f4l, Q, 0). Recently, Khastan and Nieto 111 have 
found solutions for a large enough class of boundary value 
problems with the generalized derivative. However as it is seen 
from the examples in mentioned article, these solutions are 
difficult to interpret because four different problems, obtained 
by using the generalized second derivatives, often does not 
reflect the nature of the problem. 

The second approach is based on generating the fuzzy 
solution from the crisp solution. In particular case, for the 
fuzzy initial value problem this approach can be of three 



ways. The first one uses the extension principle. In this 
way, the initial value is taken as a real constant, and the 
resulting crisp problem is solved. Then the real constant in the 
solution is replaced with the initial fuzzy value. In the final 
solution, arithmetic operations are considered to be operations 
on fuzzy numbers (fSl, ||9l). The second way, offered by 
Hiillemerier |10|, uses the concept of differential inclusion. 
In this way, by taking an alpha-cut of initial value, the given 
differential equation is converted to a differential inclusion and 
the obtained solution is accepted as the alpha-cut of the fuzzy 
solution. Misukoshi et al lITTl have proved that, under certain 
conditions, the two main ways of the approach are equivalent 
for the initial value problem. The third way is offered by 
Gasilov et al [12]. In this way the fuzzy problem is considered 
to be a set of crisp problems. 

In this study, we investigate a differential equation with 
fuzzy boundary values. We interpret the problem as a set of 
crisp problems. For linear equations, we propose a method 
based on the properties of linear transformations. We show 
that, if the solution of the corresponding crisp problem exists 
and is unique, the fuzzy problem also has unique solution. 
Moreover, we prove that if the boundary values are triangular 
fuzzy numbers, then the value of the solution is a triangular 
fuzzy number at each time. We explain the proposed method 
on examples. 

II. Fuzzy Boundary Value Problem 

Below, we use the notation u — {uL{r),uu{r)), (0 < 
r < 1), to indicate a fuzzy number in parametric form. We 
denote u = ul(0) and u = wh(0) to indicate the left and the 
right limits of u, respectively. We represent a triangular fuzzy 
number as w = {I, m, r), for which we have u — I and u = r. 

In this paper we consider a fuzzy boundary value problem 
(FBVP) with crisp linear differential equation but with fuzzy 
boundary values. For clarity we consider second order differ- 
ential equation: 

e" + ai{t)x' + a2{t)x ^ f{t) 

x{0)=A (1) 

xiT) = B 

We note that the coefficients of the differential equation are 
not necessary constant. 



Let us represent the boundary values as A — acr + a and 
B ^bcr + b, where acr and bcr are vectors with possibility 1 
and denote the crisp parts (the vertices) of A and B; a and b 
denote the uncertain parts with vertices at the origin. 

We split the problem ([TJ to following two problems: 

1) Associated crisp problem (which is non-homogeneous): 

' + ai{t)x' + a2{t)x = f{t) 

x{0) = acr (2) 

X{T) = bcr 

2) Homogeneous problem with fuzzy boundary values: 

x" + ai{t)x' + a2{t)x = 

xiO) = a (3) 

x{T) = b 

It is easy to see that, a solution of the given problem 
(H) is of the form x{t) — Xcr{t) + Xun{t) (crisp solution + 
uncertainty). Here Xcr{t) is a solution of the non-homogeneous 
crisp problem (|2]l; while Xun{t) is a solution of the homo- 
geneous problem (O with fuzzy boundary conditions. Xcr{t) 
can be computed by means of analytical or numerical methods. 
Hence, ([T]i is reduced to solving a homogeneous equation with 
fuzzy boundary conditions (|3]l. Therefore, we will investigate 
how to solve this problem. 

We assume the solution Xun of the problem (O be a fuzzy 
set X of real functions such as x{t). Each function x{t) must 
satisfy the differential equation and must have boundary values 
a and b from the sets a and b, respectively. We define the 
possibility (membership) of the function x{t) to be equal to 
the least possibility of its boundary values. 

Mathematically, the fuzzy solution set can be defined as 
follows: 

X = {x{t) I x" + ai{t)x' + a2{t)x = 0; 

a;(0) = a; x{T) = b; a ea; b eb} (4) 

with membership function 

Mx(a;(i)) = min {fia{a), fJ-^ib)} (5) 

The solution X, defined above, can be interpreted as a fuzzy 
bunch of functions. 

One can also interpret that we consider a FBVP as a set of 
crisp BVPs whose boundary values belong to the fuzzy sets a 
and b. 

A. A matrix representation of the solution in the crisp case 

Here we consider crisp BVP for second order homogeneous 
linear differential equation: 



x" + ai{t)x' + a2{t)x = 
a;(0) = a 
x{T) = b 



(6) 



Let xi{t) and X2{t) be linear independent solutions of the 
differential equation. Then the general solution is x{t) = 



ciXi{t) + C2X2{t). For Cl and C2 we have the following linear 
system 

cia;i(0) +02X2(0) = a 

cixi(r) +C2X2(r) = 6 



(7) 



Below we obtain a matrix representation for the solution of 
the BVP. We rewrite the linear system ^ in matrix form: 



M c 

a;i(0) X2(0) 
_ xiiT) X2{T) _ 
The solution of the linear system is 



where M — 



= u 








; c — 


Cl 
C2 _ 


; u = 


a 
b _ 



c = M-^ u 



(8) 



We constitute a vector-function of linear independent solutions 
s{t) = [xi{t) X2{t)]. Then the general solution can be 
rewritten in matrix form as 



Cl 
C2 



= s{t) C 



X{t) = [Xi{t) X2it)] 

Using dU we have x{t) = s(t) M^^ u, or, 

x{t) = w(t) u ~wi{t) a + W2{t) b 

where 

w(t) = s(t) M-^ 



(9) 
(10) 



B. The solution method for FBVP 

Now we show how to find X{t) (the value of the solution 
for the problem (|3]l at a time t). 

Let linear independent solutions of the crisp equation 
(l3]l, xi{t) and X2{t), be known. Then we can constitute the 
vector w (see, formula (fToli). According (|4]i and (|9]l we have: 



X = <^ x[t) = w(i) u I u : 



a ^ a\ & e 6 > (11) 



Consider a fixed time t. Put v = w(t). Then from (fTTT l we 
have: 



X(t) = i V u I u = [a 6]^ ; a e a; 6 G 6 [ 



(12) 



To determine how is the set X{t) we consider the transfor- 
mation T(u) = V u (here v is a fixed vector). One can see that 
T : i?^ ^ i?^ is a linear transformation. Therefore, X{t) is the 
image of the set B = < u = [a b] | a g a; fe G 6 > = (a, fe) 
under the linear transformation r(u). 

We remember some properties of linear transformations 

ma: 

1. A linear transformation maps the origin (zero vector) to 
the origin (zero vector). 

2. Under a Unear transformation the images of a pair of 
similar figures are also similar 

3. Under a linear transformation the images of nested figures 
are also nested. 

In addition, we shall reference a property of fuzzy number 
vectors. 

4. The fuzzy set B — (a, b) forms a fuzzy region in the 
a6-coordinate plane, vertex of which is located at the origin 



and boundary of which is a rectangle. Furthermore, the a-cuts 
of the region are rectangles nested within one another. 

The facts 1-4 allow us to derive the following conclusion. 
The vector B, components of which are the boundary values 
a and b, form a fuzzy rectangle in the afo-coordinate plane. 
The linear transformation T{u) maps this fuzzy rectangle to a 
fuzzy interval in R^. Therefore, the solution at any time forms 
a fuzzy number. 

C. Particular case when boundary values are triangular fuzzy 
numbers 

In particular, if a and b triangular fuzzy numbers, the a- 
cuts of the region B = (a, b) are nested rectangles, further- 
more, they are similar. According to the discussion above, 
their images are intervals that also are nested and similar, 
consequently, form a triangular fuzzy number X(i). Therefore, 
X{t) can be represented in the form X{t) = {x[t),0,x{t)). 
Now we investigate how to calculate x(t) and x{t). 

Let a — (a, 0,a), b — {b,0,b) and w(t) = {wi{t),W2{t)). 
Since x{t) is the maximum value among all products 
w(t)-u =awi{t) + bw2{t), we have: 

x{t) — max {awi (t) , awi (t) } + max {bw2 {t), bw2 (i) } 
x_{t) — min {awi (t) , awi (t) } + min {bw2 (t) , bw2 (t) } 

Note that an a-cut of X{t) can be determined by similarity: 

Xa{t) = \x^{t),x^{t)] = {l-a)[xit),x{t)] = {l-a)Xo{t) 

Formulas for x{t) and x{t) allow us to represent the solution 
in a new way: 

X{t) = wi(t) a + W2it)b 

where the operations assumed to be multiplication of real 
number with fuzzy one, and addition of fuzzy numbers. 

D. General case when boundary values are parametric fuzzy 
numbers 

In the general case, when a and b are arbitrary fuzzy 
numbers, the solution can be obtained by using a-cuts. Let 
Oq = [0^,0^] and ba — [60,6^. Then Ba ~ [0^,0^] x 
[ba , 5^ . By similar argumentation to the preceding case, for 
the a-cut of the solution we obtain the following formulas: 

Xa{t)^ [x^it),X^it)] 

where 

'x^{t) — max I agWi {t),'a^wi (t) ] +max \baW2{t), baW2{t)} 
Xajt) = min \ aaWi{t) ,a^wi{t)] + min {baW2{t),baW2{t)} 
Based on the formulas above we can conclude, that the 
representation for solution 



X{t) ^wi{t) a + W2{t) b 



(13) 



is valid in general. Thus, the solution, defined by formula 
dUl, becomes the same as the solution obtained from (|9]l by 
extension principle. 

Remark: The approach is valid also for the general case, 
when nth-order ?7i-point boundary value problem is consid- 
ered. 



E. Solution algorithm 

Based on the arguments above, we propose the following 
algorithm to solve the problem ([T]i: 

1 . Represent the boundary values as yl = Ucr + a and B = 
bcr + b. 

2. Find linear independent solutions xi{t) and X2{t) of 
the crisp differential equation x" + ai{t)x' + a2{t)x = 0. 
Constitute the vector-function s{t) — [xi{t) X2{t)], the 
matrix M and calculate the vector w(i) — {wi{t),W2{t)) by 
formula (fTOl i. 

3. Find the solution Xcr{t) of the non-homogeneous crisp 
problem (|2]l. 

4. The solution of the given problem ([U is 

X{t) =Xcr{t)+Wi{t)a + W2{t)b (14) 

III. Examples 
Example 1. Solve the FBVP: 

x" -3x' + 2x = At-6 

x{0) = (L5, 2, 3) (15) 

x(l) = ( 2, 3, 4) 

Solution: We represent the solution as x{t) = Xcr{t) + Xun{t). 
1 . We solve crisp non-homogeneous problem 

x" -3x' + 2x = 4:t-6 

x{0) = 2 

xll) = 3 
and find the crisp solution 

x,rit) - 2i + ^ [2(e2+* - ei+2*) + (e^* - e*)] (the thick 
line in Fig.l). 



X 3 




Fig. 1. The fuzzy solution and its a = 0.5-cut for Example 1. Dashed and 
thick bars represent the values of the fuzzy solution and its a = 0.5-cut at 
different times, respectively. 

2. We consider fuzzy homogeneous problem 

x" -3x' + 2x = 

x{0) = (-0.5, 0, 1) 

x(l) = (-l, 0, 1) 



xi{t) = e* and X2{t) 



^2t 



are linear independent solutions 



for the differential equation. Then s(t) = [e* e^'], 





■ 1 1 " 




e e 


= s 


(i) M-i 



and 

1 



o2+t 



a+2t 



o2t 



The formula (fT3] ) gives the solution of homogeneous problem: 



Xun{t) = ^((e2+* -ei+2*) (-0.5, 0, 1)+ 
(e2*_e*)(-l, 0, 1)) 



(16) 



where the arithmetic operations are considered to be fuzzy 
operations. We add this solution to the crisp solution and get 
the fuzzy solution of the given FBVP dTsl l: 



x{t) = 2t - 



(e2*-e*)(0, 1, 2)) 



) (1-5, 2, 3)+ 



(17) 



The fuzzy solution x{t) forms a band in the te-coordinate 
plane (Fig. 1). Since wi{t) > and W2{t) > for < t < T, 
the upper border of the band, x{t), becomes the solution of 
the crisp non-homogeneous problem with the upper boundary 
values A = 3 and B = A, while the lower border x_(t) 
corresponds to A= 1.5 and B_= 2: 



x{t) = 2t - 
xit) = 2t - 



1 



-((e 



2+t 



„l+2i 



)-1.5 + (e2*-e*)-0) 



e"" — e 
We can express the solution x{t) also via a-cuts, which 



are intervals Xa{t) 



at any time t. Since the 



x^{t), x^(t) 

boundary values are triangular fuzzy numbers, Xun{t) also is 
a triangular fuzzy number, say x„„(i) = (x„„(i), 0, lihm{t)). 
Consequently, an a-cut of a?„„(i) can be determined by 
similarity with coefficient (1 — a), i.e. 

Adding the crisp solution gives the a-cut of the solution x{t): 



x^{t), x^{t) 



(e^* - e*) [0, 2]) 



2+i 



a+2t 



) [1.5, 3]. 



3.5) 
1.5) 



(18) 



In Fig. 1 we show the fuzzy solution (dashed bars) and its 
a = 0.5-cut (thick bars) at different times. 

Example 2. Solve the FBVP: 

x" + 16a; = 47 - 
a;(0)-( 2, 3, 
x(2) = (0.5, 1, 

Solution: 

Associated crisp non-homogeneous problem 
mx = 47 - 8^2 
x{0) = 3 
x(2) = 1 
has the solution Xcr{t) == 3 — 0.5t^ (thick line in Fig. 2). 

To find the uncertain part of the fuzzy solution, Xun{t), we 
solve fuzzy homogeneous problem 



X 2 




Fig. 2. The fuzzy solution and its a = 0.6-cut for Example 2. Dashed and 
thick bars represent the values of the fuzzy solution and its a = 0.6-cut at 
different times, respectively. 



x" + 16x = 
x{0) = ( -1, 0, 0.5) 
x{2) = (-0.5, 0, 0.5) 
a:;i(t) = cos4i and X2{t) = sin4i are linear independent 
solutions for the differential equation. Then 

' 1 



s(t) = [cos4i sin4i], M 



and 



cos 8 sin 8 
w = s(t) A/-1 = 5jig [sin(8 -At) sin4t]. 

Using the formula (fT3] l we obtain the solution of homo- 
geneous problem and adding the crisp solution we get the 
solution of the given FBVP (fTsT l: 

x{t) = 3 - 0.5i2 + 1 (sin(8 - At) (-1, 0, 0.5)+ 



sin4i (-0.5, 0, 0.5)) 



(19) 



Fuzzy solution generates a band in te-plane (Fig. 2). Unlike 
Example 1, the functions wi{t) and ■W2(t) takes both positive 
and negative values in the interval < i < T. Because of 
that, in generation of upper and lower borders of the band, a 
and a, b and b take charge in alternately. 



IV. Conclusion 

In this paper we have investigated the fuzzy boundary value 
problem as a set of crisp problems. We have proposed a 
solution method based on the properties of linear transfor- 
mations. For clarity we have explained the proposed method 
for second order linear differential equation. We have shown 
that the fuzzy solution by our method coincides with the 
solution by extension principle. We are planning to make a 
comparative analysis between the proposed method and the 
method including generalized Hukuhara derivative in future. 
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